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.. , ^ c#nr i- xhis tvDe of formulation provides both general insight into PZT 

\ nonlinear lumped-parameter model ot a piezoelectric ■ sta JP „ as a specl f ic causa! mathematical representation 

for microrobotic applications requiring accurate position and/or 
force control. In formulating this model, the authors propose a 
generalized Maxwell resistive capacitor as a lumped-parameter 
causal representation of rate-independent hysteresis. Model for- 
mulation is validated by comparing results of numerical simu a- 
tions to experimental data. Validation is followed by a discussion 
of model implications for purposes of actuator control. 


behavior as well as a specific causai mauicuuuw. ■ 

for purposes of model-based control system analysis and design. 

The Anatomy of a Piezoelectric Stack Actuator 

Piezoelectric ceramics transduce energy between the electri- 
cal and mechanical domains. Application of an electric field 
across the ceramic creates a mechanical strain, and in a simi ar 
manner, application of a mechanical stress on the ceramic in- 
duces an electrical charge. The fundamental component of a pie- 
zoelectric stack actuator is a wafer of piezoelectric maten 
sandwiched between two electrodes. Prior to fabrication, the wa- 
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Micromanipulator TOlnua-j JLo n significant pieroelccmc effect in this direction only A typical 
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rally interactive micromanipulation, require control actuators 
that can provide both accurate position tracking performance in 
addition to robustly stable force control [2, 7, 13]. These objec- 
tives are significantly compromised by the presence of backlash 
and Coulomb friction in the control plant [4, 14, 16], the effects 
of which are exaggerated at small scales [5]. Since piezoelectric 
stack actuators are monolithic and have no sliding or rolling 
parts, they exhibit no significant mechanical sticuon or backlash. 
Additionally, a typical lead-zirconate-titanate (PZT) piezoelec- 
tric stack actuator can perform step movements with a resolution 
on the order of a nanometer, has a bandwidth on the order of a ki- 
lohertz, and can provide mechanical power on the order of sev- 
eral watts. Consequently, piezoelectric ceramics are well suited 
. . „ c — nevices. 


ter IS poiiuiicu uiHOAifti m^.A & 

significant piezoelectric effect in this direction only. A typical 
piezoelectric stack actuator is formed by assembling several of 
the wafer elements in series mechanically and connecting the 
electrodes so that the wafers are in parallel electrically, as ll us- 
trated in Fig. 1 . The nominal quasi-static behavior of a piezoelec- 
tric stack actuator is a steady-state output displacement that is 
monotonically related to the voltage input. ... _ 

Existing Descriptions of Piezoelectric Behavior ~ 

The most widely recognized description of piezoelectric ce- 
ramic behavior was published by a standards committee of the 
IEEE Ultrasonics. Ferroelectrics, and Frequency Control Soci- 
ety originally in 1966 and most recently reused m 1987 [15J. 
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and due to the absence of Coulomb friction and backlash, are es- 
pecially well suited for use in dynamically interactive microma- 

n ' P The use ofl PZT stack actuators for accurate and stable control 
of manipulator position and/or force is greatly facilitated by 
model-based control system analysis and design. The purpose ot 
the model presented herein is to map the relationship between 
voltage and charge at the electrical port of the PZT to force and 
displacement at the mechanical port in a lumped parameter form 
that can be represented by a set of ordinary differential equations. 


ims conuuiucc . . 

scribing piezoelectric continua which form the basis for the 
model of piezoelectric behavior that is presently in general use. 
The linearized constitutive relations are typically represented in 
a compressed matrix notation as follows. > , 

(1) 


D, = dj, + t T ,A 


( 2 ) 
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where S represents the strain tensor. / is the elastic compliance 
matrix when subjected to a constant electrical field. T represents 
the stress tensor, d is a matrix of piezoelectric material constants, 
E is the electric field vector. D the electric displacement vector, 
and e r the permittivity measured at a constant stress. The com- 
pressed notation eliminates redundant terms by representing the 
symmetric stress and strain tensors with single column vectors 
that incorporate elements representing both the diagonal and 
off-diagonal tensor terms. These equations essentially state tha 
the material strain and electrical displacement ( ch “ r 8 e P c 'j 
area) exhibited by a piezoelectric ceramic arc both Iincarty af^- 
fectcd by the mechanical stress and electrical held to which the 
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Fig. I . Illustration of a piezoelectric stack actuator. 


ceramic is subjected. Aside from the awkward notation, these 
linearized constitutive relations fail to explicitly describe the 
nonlinearities that are present in all piezoelectric ceramics. Ad- 
ditionally, the derivation assumes a purely conservative energy 
field, and thus fails to describe the dissipative behavioral aspects 
of the ceramic. 

Several authors have utilized these linearized constitutive re- 
lations to derive a piezoelectric actuator model [6, 10, 12]. Given 
the aforementioned limitations of the IEEE model, however, the 
resulting descriptions are too approximate for use in developing 
accurate position or force control. Leigh and Zimmerman took 
an iterative algorithmic approach to describing the nonlinear be- 
havior of piezoelectric ceramics [11], Though capable of de- 
scribing the hysteretic behavior of piezoelectrics, this approach 


is neither real-time nor energy-based, and thus provides little 
physical insight into actuator behavior and has limited use with 
respect to model-based control system design. Jouaneh and Tian 
incorporated a time-dependent nonlinear state variable, a 
method originally developed by Bouc for a structures applica- 
tion, to describe the nonlinear voltage/displacement behavior of 
piezoelectric actuators [1,8]. The model of Jouaneh and Tian is 
nicely posed, but as an electromechanical model is incomplete, 
since it represents both the electrical and mechanical ports of the 
PZT as signal ports (i.e., the model implicitly assumes no dy- 
namic interaction with the environment). Utilizing piezoelectric 
ceramic as an actuator explicitly requires that it perform work on 
its environment. An accurate behavioral description must there- 
fore characterize the input and output ports of the PZT as power 
ports, and must therefore describe relationships involving me- 
chanical force and electrical charge as well as electrical voltage 
and mechanical displacement. 

The model presented herein offers a complete behavioral de- 
scription that is based neither on the IEEE constitutive relations, 
nor on algorithm, but rather is formulated by describing the ob- 
served static and dynamic behavior of a piezoelectric actuator 
with a lumped-parameter energy-based representation. 

Model Formulation 

Since piezoceramic is a known dielectric, one would expect a 
PZT stack actuator to exhibit capacitive behavior. The electrical 
behavior of the actuator, however, is significantly more complex. 
For purposes of controller design, one of the most inconvenient 
aspects of the actuator behavior is the rate-independent hystere- 
sis exhibited between voltage and displacement as well as be- 
tween force and displacement, as shown in Fig. 2. If not 
specifically addressed, this type of behavior can cause closed- 
loop limit cycling and possibly instability. Experimental obser- 
vation indicates that the rate-independent hysteresis exhibited in 
Fig. 2 is not present between the endpoint displacement of the 



Fig. 2. Measured quasi-static relationships between applied voltage and endpoint displacement and between applied force and endpoint 
displacement. 
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f/c 3 Column (a)- Measured charge input versus endpoint displacement output for PZT stack actuator for 100 and 1000 

inputs, tolumn (b . ): Frequency response of endpoint displacement output to charge input. The marks represent measured data points and the 

continuous lines represent a linear second-order system. 


PZT stack actuator and the net electrical charge delivered to the 
actuator. Additionally, dynamic observation indicates that end- 
point displacement as a function of electrical charge is well ap- 
proximated by second-order linear dynamics, as shown in the 
measured data of Fig. 3. 

The quasi-static force-displacement relationship of a PZT 
stack actuator is shown in Fig. 4. As illustrated in the figure, 
rate-independent hysteresis is observed only when the electrode 


leads are shorted. When the leads are open and current cannot 
flow through the ceramic, the actuator exhibits no static hystere- 
sis. This evidence, along with the absence of static hysteresis be- 
tween displacement and charge illustrated by Fig. 3(a). suggests 
that the rate-independent hysteresis lies solely in the electrical 
domain between the applied actuator voltage and resulting 
charge. 
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PZT Force/Displacement Behavior 



Fig. 4. Quasi-static force displacement relationship fora PZT stack 
actuator with open leads ( dashed line) and with shorted leads (solid 
line). 


Modeling Hysteresis 

Generalized hysteretic behavior is characterized by two vari- 
ables related by a non-single-valued function, as illustrated by 
the voltage/displacement relationship in Fig. 2. Examples of 
such behavior include the relationship between the stress and 
strain in the elastic -plastic deformation of a material, and the re- 
lationship between magnetic field strength and flux density in a 
hard magnetic material. In these and in most physical systems, 
this hysteretic behavior is the result of energy storage that is fun- 
damentally coupled to rate-independent dissipation. Mechani- 
cally, such behavior could be modeled by the combination of an 
ideal spring, which represents pure energy-storage, coupled to a 
pure Coulomb friction element, which represents rate- 
independent dissipation. This analogy is the basis for describing 
the static hysteresis exhibited by the piezoelectric actuator. Fig. 5 
illustrates a single elasto-slide element which consists of a mass- 
less linear spring and a massless block that is subjected to Cou- 
lomb friction. The constitutive behavior of this element can be 
described by: 


[ i{x-x h ) \f\k(x-x h ]< f 

F = i / 

/ sgn(.t) and x b - x - — sgn(i) else 
l k (3) 

where .r is the input displacement, F is the output force, k is the 
stiffness of the spring,/ = fJjV is the breakaway friction force of 
the block, and xh is the position of the block. For a displacement 
input of sufficient amplitude, the relationship between the ap- 
plied force and the endpoint displacement will exhibit rudimen- 
tary hysteretic behavior. If several of these elasto-slide elements 
are put in parallel, each subjected to an incrementally larger nor- 
mal force (i.e., each has an incrementally larger breakaway 
force), the simple relationship of Fig. 5 becomes a piecewise lin- 
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Fig. 5. Force-displacement behavior of a single elasto-slide element 
consisting of a massless linear spring and a massless block 
subjected to Coulomb friction. 


ear approximation of rate-independent hysteresis, as illustrated 
in Fig. 6. This construction was initially formulated by the 
mathematician and physicist James C. Maxwell in the mid : 
1 800s, and in the limit as the number of elasto-slide elements be- 
comes infinite, the model is referred to as Generalized Maxwell 
Slip [9], The constitutive behavior of the system in Fig. 6 con- 
taining n elasto-slide elements can be described by: | 


*.(- l - J O if K( r "^)| < / 

f, sgn(.t) and x b ~ x ~~ sgn (i) else 


(4) 


F-iF. 


(5) 


where .r is the input displacement, F is the output force, and Ft, ku 
fi , and.r^ are the output force, spring stiffness, breakaway force, 
and block position, respectively, of the I th elasto-slide element. 

The Maxwell model of rate-independent hysteresis can be pa- 
rameterized directly from a piecewise linear fit of the rising 
curve of the hysteresis from a relaxed state, as illustrated in Fig, 
7. Note that complete parameterization of the entire hysteretic 
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Fig. 6. Force-displacement behavior of several elasto-slide elements 
in parallel, each subjected to increasing normal forces. 

behavior requires fitting only the initial rising curve. A rising 
curve fit by n linear segments will require n elasto-slide ele- 
ments, and will thus require In parameters. For a curve fit by n 
linear segments, the slope of the j lh segment (the localized stiff- 
ness) is given by: 

w (6) 

where to represents the stiffness of the i th spring in the Maxwell 
model. This can be expressed in a matrix form as: 

s = A k (?) 

where s is a at x l vector of the segment slopes, A is an at x ai upper 
triangular matrix of ones, and k is an ai x l vector of Maxwell 


Fig. 8. Maxwell model parameters corresponding to hysteresis and 
piecewise linear fit of Fig. 7. 

model spring stiffness, as illustrated in Fig. 8. Rearranging gives 
an explicit expression for the spring stiffnesses: ,. 

k » A*‘s. _<8) 

•*- r 

The location of each respective segment, defined by Xj in the five 
segment fit of Fig. 7, can be expressed as: /jfc. 
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Fig. 9. Schematic representation of the PZT stack actuator model \ 
The capacitor and resistor contained within the MRC element are 
both nonlinear elements . 


where fj is the breakaway friction force of the j 01 block in the 
Maxwell model. Rearranging gives an explicit expression for the 
breakaway forces: 


f = K x (10) 

where f is an n x 1 vector of the breakaway forces, K is an n x n di- 
agonal matrix of the spring stiffnesses, and x is an n x 1 vector of 
the segment locations as shown in Fig. 7. Note that since the 
Maxwell model is a piecewise linear approximation, increasing 
the number of blocks will increase the model accuracy. The 
number of elasto-slide elements, however, has no bearing on the 
order of the model, since the blocks are all massless and the 
springs are in parallel. 

Though ostensibly a mechanical formulation, the energy- 
based constitutive relations of the Maxwell slip model are not 
domain-specific, and can therefore represent any rate- 
independent hysteretic relationship between a generalized force 
and generalized displacement in a lumped parameter causal 
form. Consequently, in addition to force and displacement, the 
generalized Maxwell model can represent rate-independent hys- 
teresis between voltage and charge, temperature and entropy, 
and magnetomotive force and magnetic flux. 

The PZT stack actuator model resulting from the aforemen- 
tioned observations is shown in schematic form in Fig. 9. The 
generalized Maxwell resistive capacitance, which is represented 
by the MRC element, resides in the electrical domain and there- 
fore relates the element’s electrical voltage to charge. The PZT 
model has two ports of interaction, a voltage -current port on the 
electrical side and a force-velocity port on the mechanical side. 
With respect to the mechanical side of the transformer, since the 
actuator model is concerned only with endpoint displacement in 
a frequency band within the first mechanical mode of vibration, 
the piezoelectric stack is assumed to have a lumped mass and a 
linear material stiffness and damping. The behavior of the actua- 
tor is therefore described by: 

q = Tx+Cv t (11) 

= K + K ( 12 ) 


F t = 7V, (14) 

mx + bx + kx = F t + f rU (15) 

where q is the total charge in the ceramic, T is the electrome- 
chanical transformer ratio, x is the stack endpoint displacement, 
C is the linear capacitance in parallel with the transformer, v, is 
the back-emf from the mechanical side, v m is the actuator input 
voltage, \?rc is the voltage across the Maxwell capacitor (which as 
indicated by Equation ( 1 3) is a function of q ), F t is the transduced 
force from the electrical domain, m, b , and k are the mass, damp- 
ing, and stiffness of the ceramic, and is the force imposed 
from the external mechanical load. The function of equation ( 1 3) 
relates the voltage across the MRC element to the charge in the 
ceramic as follows: 


Table 1. Model Parameters Utilized in Simulation of the NEC 
Model #AE0505D16 Piezoelectric Actuator 

Model Parameter 

Symbol 

Numerical Value 

Mass 

m 

0.00375 kg 

Stiffness 

k 

6 x 10 6 N/m 

Damping 

b 

150 N-s/m 

Linear Capacitance 

C 

1.2 x 10" 6 F 

Transformer Ratio 

T 

1 10 am 


V=. 


(*-*>): 




if 




< v,. 


v, sgn(j') and q b =q- C,v, sgn(<) else ( j 



*1 




i 



where q is the charge in the ceramic, v r c is the voltage across the 
Maxwell capacitor, and all other variables are the electrical ana- 
logs of those given in Equations (4) and (5). 


Simulation and Experiment 
The model was parameterized for a commercially available 
piezoelectric stack actuator (NEC model #AE0505D16). This 
actuator operates at input voltages between zero and 150 volts (in 
the direction of wafer polarization), which corresponds to an 
endpoint displacement range of approximately 20 microns. 
Upon measuring the mass m of the PZT, the mechanical stiffness 
k and damping b can be determined by observation of the 
charge-displacement dynamics that are shown in Fig. 3. Meas- 
urement of the open-lead stiffness and knowledge of the DC gain 
between charge and displacement are sufficient to determine the 
linear electrical capacitance C and the transformer ratio T The 
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Equation (8) is replaced with the 
inverse of a capacitance and each 
breakaway force of Equation (10) 
is replaced with an analogous 
saturation volume: 


— = * 
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V =/■ 


(18) 
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Fig. 10. Measured (solid line) and simulated ( dashed line) endpoint displacement response of the 
PZT to a 90 volt 100 Hz triangle-wave input. The difference between the measured and simulated 
data is difficult to discern. 


model parameters utilized for the simulations that follow are 
given in Table 1 . 

The Maxwell capacitor parameters, as given in Table 2, were 
determined by piecewise linear fit of the measured voltage ver- 
sus charge data. More specifically, the electrical parameters of 
Equations (16) and (17), C\ and vi , were found using the meth- 
ods of Equations (6) through (10), where each stiffness of 
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The parameters defining the 
Maxwell capacitor for the simula- 
tions presented herein are given in 
Table 2. As indicated by the table, 
these simulations were run with ten 
elasto-slide elements in the gener- 
alized Maxwell resistive capacitor. 
As previously stated, the Maxwell 
model is essentially a piecewise lin- 
ear approximation of rate- 
independent hysteresis where the 
number of elasto-slide elements 
corresponds to the number of linear 
segments. Increasing the number of 
blocks therefore increases the 
model accuracy. The number of 
blocks required to accurately 
model a given hysteresis thus de- 
pends on the curvature of the hys- 
teresis and the required degree of 
model accuracy. A real-time appli- 
cation may be better served by 
fewer elements; since the number 
of elements does not affect the or- 
der of the model, however, the in- 
creased computational overhead 
from added elements is minimal..*. 

Fig. 10 shows the measured and 
simulated endpoint displacement 
response of the modeled actuator to 
a 30, 60, and 90 volt 100 Hz 
triangle-wave voltage input in the 
absence of an external mechanical 
load. As indicated by the plots, the 
model faithfully represents the 
voltage input to endpoint displace- 
ment output behavior of the PZT 
actuator. Note in the plots that the 
displacement output evinces distor- 
tion on both the rising and falling 
slopes and maintains an 
amplitude-dependent offset. Fur- 
ther observation indicates that this is not typical dynamic dis- 
tortion. There is no discernible phase lag between the input and 
output, as indicated by the relative positions of the waveform 
peaks, and there are no significant filtering effects, evidenced 
by the fact that the peaks are neither rounded iu>r otherwise dis- 
torted. These observations indicate the existence ot a rate- 
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Tiible 2. Generalized Maxwell Capacitor 
Parameters Utilized in the Model* 

Element 

Stiffness 

Break Force 

i 

ki (I/O 

fi (Vf) 

1 

2.0 

0.2 

2 

0.6 

0.3 

3 

0.3 

0.3 

4 

0.26 

2.6 

5 

0.06 

0.9 

6 

0.1 

2.0 

7 

0.05 

1.5 

8 

0,03 

1.2 

9 

0.1 

7.0 

10 

0.5 

80.0 

* The stiffnesses and breakaway forces correspond with the me- 
chanical schematic of Fig. 6. 


independent hysteresis. This hysteretic behavior is clearly dis- 
played in Fig. 1 1, which shows both the measured and simulated 
data of Fig. 10 plotted as voltage input versus displacement out- 
put. Any dynamic system will exhibit an ellipsoid-shaped input- 
output hysteresis, provided the input is of sufficient frequency to 
create discernible phase lag in the output. This type of dynamic 
hysteresis, characterized by a smooth curve relating the output to 
input, is a result of phase lag between the input and output (such 
as that shown in Fig. 3). The rate-independent hysteresis of Fig. 
10, however, is distinguished from a dynamic hysteresis by the 
distinct discontinuities exhibited at both extremes. Model accu- 
racy is additionally demonstrated by the data of Fig. 12, which 
shows the measured and simulated response of the actuator to a 
linearly decaying 100-hertz sinusoidal voltage input. r 

Fig. 13 shows the simulated displacement response to a 10 Hz 
sinusoidal external force for cases with the electrical leads open 
and shorted. The simulated behavior accurately reflects the qua- 
sistatic measured behavior illustrated by Fig. 4. In addition to pro- 
viding a structure for numerical simulation, the lumped-parameter 
model also provides insight into how the electrical properties of 
the PZT reflect into the mechanical domain. As shown by the 
measured and simulated data of Figs. 4 and 13, respectively, the 
PZT actuator exhibits significantly greater stiffness when the 
leads are open than when the leads are shorted. This change in 
stiffness can be demonstrated by observing that both the linear and 
the Maxwell capacitors in the electrical domain reflect as stiff- 
nesses in the mechanical domain. This is similar to the electrical 
resistance of a DC motor appearing as mechanical damping when 
the motor leads are shorted. Linearizing the Maxwell capacitor 
and deriving expressions for mechanical stiffness in both the open 
and shortcd-lead configurations yields 

T 2 

k t = k + — 

C (20) 


I 
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Hysteresis for 60 Volt Input 
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Hysteresis for 90 Volt Input 



Fig, II. Measured ( solid line) and simulated (dashed line) voltage 
versus displacement for a 90 volt 100 Hz triangle-wave voltage 
input. The difference between the measured and simulated data is 
difficult to discern. 
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Simulated Force/Displacement Behavior 



Displacement (Microns) 

Fig. 13. Simulated PZT behavior showing force versus endpoint 
displacement for a 10 Hz sinusoidal external force input for cases 
with the electrical leads open ( dashed line) and shorted ( solid line). 



Fig. 14.1EEE lumped parameter circuit model of PZT near isolated 
resonance . 




T 2 


C.+C 


( 21 ) 


where fc# is the open lead stiffness, k , is the shorted lead stiffness, 
k is the mechanical domain stiffness, C is the linear electrical ca- 
pacitance, T is the transformer ratio, and C„, is the linearized 
Maxwell capacitor. The inclusion of Cm in the second term of the 
latter expression clearly indicates the decreased shorted lead 
stiffness. 


Comments on the IEEE Resonance Model 

In addition to formulating the previously described linearized 
constitutive relations for piezoelectric continua, the IEEE stan- 
dards publication [15] also offers a circuit representation of a 
PZT near an isolated resonance, shown in Fig. 14. The standards 
publication does not relate the resonance circuit to the constitu- 
tive relations and does not offer any explanation for the presence 
of inductive or resistive effects in what is, according to the con- 
stitutive equations, a multi-port capacitive medium. The model 
derived herein elucidates the structure of the resonance model 
and accounts for the inductive and resistive elements. The induc- 
tor, resistor, and capacitor (L\ t i?i, and C i, respectively) in the 
IEEE resonance circuit of Fig. 14 are not fundamentally electri- 
cal behaviors, but rather are the apparent behavior of the first 
mode of mechanical vibration, as seen through the electrical 
port. Reflecting the mechanical elements of the lumped parame- 
ter model to the electrical side gives circuit element values of 

(23) 


(24) 


(25) 

(26) 

This would indicate that driving a typical mass-spring- 
damper external load would effectively increase L\ and R\ y de- 
crease Cu and leave Co unaffected. . 

Implications for Actuator Control 
Several analogies can be drawn between control of a DC per- 
manent magnet motor and control of a PZT stack actuator. Fust, 
controlling charge into a PZT actuator is in many ways analo- 
gous to controlling current into a DC motor. The relationship be- 
tween input current and output displacement of a DC motor can 
be characterized by linear second-order dynamics (one pole 
serving as a pure integrator). In a similar manner, the relationship 
between input charge and output displacement of a PZT stack 
can also be characterized by linear second-order dynamics. Sec- 
ond, just as the output shaft toque is direedy proportional to in- 
put current in a DC motor when the shaft is locked, the output 
force of a PZT stack is similarly proportional to input charge 
when the stack is locked. Continuing this analogy, controlling a 
DC motor with voltage as input introduces the first-order dynam- 
ics of the stator coil inductance and resistance into the system. In 
a similar manner, controlling a PZT stack with voltage as input 
introduces the first-order dynamics of the Maxwell resistive ca- 
pacitor into the system, which unlike the inductor/resistor dy- 
namics of the DC motor, are nonlinear. 

Summary and Conclusions 
The model presented accurately represents the behavior of a 
piezoelectric stack actuator in a lumped-parameter real-time rep- 
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resentation, and can therefore be utilized for purposes of model- I l () l l.cc. T.J. Moon, and G.Y. Masada, "Extended Bond Graph Reticula- 

based control analysis and design. The static hysteresis evinced ,inn 1,1 1'icznclectric Continua.” ASME Journal of Dynamic Systems, Meets - 

by the PZT actuator was identified as energy storage coupled to umnent, tmti Control, vol, 1 17, no. 1. March, pp. 1-7, 1995. 

rate-independent dissipation, and was faithfully represented by a 1 1 1 1 T I). Leigh and D.C. Zimmerman. “An Implicit Method for the Nonlin- 

generalized elasto-slip model. Despite the presence of this non- car Modelling and Simulation of Piezoccraniic Actuators Displaying Hys- 

linearity, the relationship between charge delivered to the PZT icresis.” ASME Smart Structures and Materials. AD- vol. 24. pp. 57-63, 

and endpoint displacement of the stack was observed to have 1 1,1,1 

simple second-order linear characteristics. Design and analysis 1 1 2 1 C. Liang, F. Sun. and C.A. Rogers. “Coupled Electro-Mechanical 

of closed-loop actuator control can therefore be achieved either Analysis of Piezoelectric Ceramic Actuator- Driven Systems: Determination 

by incorporating nonlinear techniques with actuator voltage as of the Actuator Power Consumption and System Energy Transfer.** Proceed - 

command, or by utilizing linear techniques with actuator charge "w of the SPIE Smart Structures and Intelligent Systems . vol. 1917, pp. 

as command. 2X6-298. 1993. 
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